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Abstract 

We make a new multivariate generalization of the type A monomial space of a 
single variable. It is different from the previously introduced type A space of several 
variables which is an sl(M + 1) module, and we thus call it type A'. We construct 
the most general quasi-solvable operator of (at most) second-order which preserves 
the type A' space. Investigating directly the condition under which the type A' 
operators can be transformed to Schrodinger operators, we obtain the complete list 
of the type A' quasi-solvable quantum many-body systems. In particular, we find 
new quasi-solvable models of deformed Calogero-Sutherland type which are different 
from the Inozemtsev systems. We also examine a new multivariate generalization of 
the type C monomial space based on the type A' scheme. 
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Calogero-Sutherland models 
PACS: 03.65.Ge, 02.30. Jr 



1 Introduction 



It is widely known that most of the quasi-exactly solvable quantum one-body 
Hamiltonians, for which we can obtain a part of the exact eigenvalues and 
eigenfunctions in closed form [1,2], have the underlying sl(2) structure first 
discovered in Ref. [3]. Recently, it was found in Ref. [4] that the well-known 
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exactly solvable M-body Calogero-Sutherland models [5,6] have a similar Lie- 
algebraic structure of sl(M + 1). After the several new discoveries of quasi- 
exactly solvable Hamiltonians having the same sl(M + 1) Lie-algebraic struc- 
ture [7,8,9,10], the complete list of the quantum many-body systems which 
admit the sl(M + 1) structure was obtained in Ref. [11]. It turns out that all 
of them are of Inozemtsev type, which were originally found to be classically 
integrable [12,13,14]. The common feature of all the models which have an un- 
derlying Lie-algebraic structure is that they leave a finite dimensional module 
of the Lie algebra invariant. 

On the other hand, Post and Turbiner studied a classification of linear differ- 
ential operators of a single variable which have a finite dimensional invariant 
subspace spanned by monomials [15]. According to the results in Ref. [15] 
combined with the discussion in Ref. [16], there are essentially three different 
spaces of monomial type preserved by second-order linear differential opera- 
tors, except for a few special cases. Later, they were dubbed type A, B, and C, 
respectively, according to the corresponding types of A/"-fold supersymmetry 
[17]. The type A space corresponds to the sl(2) module investigated in Ref. [3] 
and second-order linear differential operators preserving it are expressed as 
quadratic forms of the 5 1(2) generators represented by first-order linear differ- 
ential operators. The sl(M + 1) module preserved by the quantum Inozemtsev 
systems was then regarded as a natural generalization of the type A monomial 
space of a single variable to several variables [18]. 

However, the other two spaces, type B and C, are not Lie-algebraic modules 
and linear differential operators preserving the type B or C spaces are not 
given through the universal enveloping algebra of any Lie algebras. Recently, 
we have successfully generalized the type C monomial space of a single vari- 
able to several variables and constructed the most general second-order linear 
differential operator preserving it [18]. Throughout our experiences in study- 
ing the latter problems and searching for a natural generalization of the type 
B monomial space of a single variable to several variables (cf. Section 8 in 
Ref. [18]), we are convinced that there are much more varieties of monomial 
type space of several variables than those of a single variable which can be pre- 
served by linear differential operators if we do not restrict ourselves to study 
such an operator that admits an underlying Lie-algebraic structure. 

In this article, we show that by making another generalization of the type 
A monomial space of a single variable, which is different from the sl(M + 
1) module investigated fully in Ref. [11], we obtain a new family of quasi- 
solvable quantum many-body systems of deformed Calogero-Sutherland type. 
The obtained models turn to have in general M-body interaction terms and 
are thus different from the Inozemtsev systems. 

The article is organized as follows. In the next section, we summarize some 
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definitions related to the concept of quasi-solvability in order to avoid ambi- 
guity. In Section 3, we introduce a new generalization of the type A monomial 
space of a single variable to several variables, which we shall call type A', 
and briefly discuss some important property of the type A' space such as the 
invariance under linear transformations. In Section 4, we construct the most 
general (at most) second-order linear differential operator which preserves the 
type A' space. In Section 5, we examine the condition under which the general 
type A' quasi-solvable operators can be transformed to Schrodinger operators. 
Utilizing the invariance under the linear transformations, we fully classify the 
possible quantum Hamiltonians preserving the type A' space in Section 6. In 
Section 7, we make a new multivariate generalization of the type C monomial 
space of a single variable, which we shall call type C, based on the type A' 
space. We then construct the most general type C quasi-solvable operator of 
(at most) second-order as well as the type C gauged Hamiltonian. Finally in 
Section 8, we summarize and discuss the obtained results in combination with 
the results in the type A and C cases. Some useful formulas are summarized 
in Appendix A. 



2 Definition 

First of all, we shall give the definition of quasi-solvability and some notions 
of its special cases based on Refs. [11,19]. A linear differential operator H of 
several variables q = • • • , ?m) is said to be quasi-solvable if it preserves a 
finite dimensional functional space V/y whose basis admits an analytic expres- 
sion 4>i(q) in closed form 2 : 



HV^cVx, dim Vat = n(Af) < oo, Vx = (Ml), ■ ■ • > 0n(jv)(g))- (2.1) 



An immediate consequence of the above definition of quasi-solvability is that, 
since we can calculate finite dimensional matrix elements Sk,i defined by, 



we can diagonalize the operator H and obtain its spectra in the space VV, at 
least, algebraically. Furthermore, if the space Vat is a subspace of a Hilbert 
space L 2 (S) (S C M M ) on which the operator H is naturally defined, the 
calculable spectra and the corresponding vectors in VV give the exact eigen- 
values and eigenfunctions of H, respectively. In this case, the operator H is 

2 The latter restriction has been sometimes missed in the literature. Without it, 
however, arbitrary linear operators would be quasi-solvable unless their spectrum is 
empty. 



n(Af) 

H(j) k = S >*,i4>i 




(2.2) 



i=i 
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said to be quasi-exactly solvable (on S) 3 . Otherwise, the calculable spectra 
and the corresponding vectors in Vjv only give local solutions of the charac- 
teristic equation of H . This important difference has been sometimes missed 
in the literature. 

A quasi-solvable operator H of several variables is said to be solvable if it 
preserves an infinite flag of finite dimensional functional spaces VV, 

Vi C V 2 C • • • C Vat C • • • , (2.3) 

whose bases admit analytic expressions in closed form, that is, 

HV Af cV M , dim Vat = n (AO < oo, V M = (<f>i (?),•••, <l>n(N)(q)), (2.4) 

for Af = 1, 2, 3, . . .. Furthermore, if the sequence of the spaces (2.3) defined on 
Scl M satisfies, 

V^(S)^L 2 (S) (jV-00), (2.5) 
the operator H is said to be exactly solvable (on S). 



3 A New Generalization of the Type A Monomial Space 

In this article, we shall consider a quantum mechanical system of M identical 
particles on a line. The Hamiltonian is thus given by 

1 M d 2 

^-2^ + %--?m), (3.1) 

where the potential has permutation symmetry: 

V(...,q i ,...,q j ,...)=V(...,q j ,...,q i ,...) Vi ^ j . (3.2) 

To construct a quasi-solvable operator of the form (3.1), we shall follow the 
three steps after Ref. [11], namely, i) a gauge transformation on the Hamilto- 
nian (3.1): 

H = e w ®He- w ( q \ (3.3) 

3 A domain S is not necessarily a subspace of the real space M M if the operator 
under consideration is non-Hermitian. Indeed, a couple of quasi-solvable one-body 
Hamiltonians which are not quasi-exactly solvable on any subspaces of the real space 
M. are shown to be quasi-exactly solvable on the subspaces of the complex space C 
incorporating with the 'PT-symmetric boundary conditions [20,21]. 
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ii) a change of variables from qi to z\ by a function z of a single variable: 

Zj = (3.4) 

and iii) the introduction of the elementary symmetric polynomials of z^ defined 
by 

M 

°k{z) = J2 z h--- z ik (k = l,...,M), (7 = 1. (3.5) 

h<—<ik 

Due to the permutation symmetry of the original Hamiltonian (3.1), the 
gauged Hamiltonian (3.3) can be completely expressed in terms of the ele- 
mentary symmetric polynomials (3.5). In this article, a second-order linear 
differential operator is called a gauged Hamiltonian if it can be transformed to 
a Schrodinger operator by means of a combination of gauge transformations 
and change of variables. 

The next task is to choose a vector space to be preserved by the gauged 
Hamiltonian (3.3). The type A monomial space of a single variable z is defined 
by [17] 

Vj^l,*,...,*"- 1 ). (3.6) 

It is an sl(2) module and the foundation of the sl(2) construction of quasi- 
solvable models in Ref. [3]. In our previous paper [18], we identified the fol- 
lowing space as a generalization of the type A space of a single variable to 
several variables: 

M . 

m G Z> , ]Tn* < A/"- l). (3.7) 
i=i ' 

It is indeed a natural generalization since it provides an sl(M + 1) module for 
each M and is the foundation of the sl(M + 1) construction of quasi-solvable 
quantum many-body systems in Ref. [11]. However, the spaces (3.6) and (3.7) 
have a different character, that is, the elements of the latter space of multi- 
variable are not necessarily polynomials of degree less than M in the variables 
Zi, in contrast to the former space of a single variable. Hence, another natural 
generalization would be such that any element of a generalized space is a 
polynomial of degree less than M in Zi. Each elementary symmetric polynomial 
<7k is a polynomial of degree k in Zi and thus the latter generalization can be 
realized by the following vector space: 

m . 
G Z> , ^in, <N- l). (3.8) 
i=i ' 

Obviously, the space (3.8) also reduces to the single-variable type A space 
(3.6) when M — 1, but is different from Eq. (3.7). We thus call the space 



Kaam 



n M 
°M 



V N;M — { a l ■ ■ ■ °. 
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(3.8) type A'. In this article, we investigate linear differential operators which 
preserve the type A' space. 



In contrast to the fact that the type A space (3.7) is invariant under the 
GL(2,M) linear fractional transformation [11]: 

M 

V { m\z] „ n(7* + Sf-^jPW = V#>[*], (3-9) 
i=i 

induced by 

ry Z . _|_ R 

Zi^Zi= -^—^ (a, /3,7,5ef;A = a5-/57^0), (3.10) 

the type A' space (3.8) does not have the full GL(2, M) invariance for M > 1. 
It is easily read from the transformation of the elementary symmetric polyno- 
mials (3.5) under the special projective transformation Z{ \— > Zi — 1/zf. 

a k (z) h-> (Tk{z~ l ) = a M ^ k ( y z)a M (zy 1 . (3.11) 

Hence, the special projective transformation interchanges the role of the vari- 
ables a k and o M _ k in the space (3.8) and cannot keep the condition Y^iLi i> n i < 
N — 1 unchanged. However, the other elements of the GL(2,M) transforma- 
tions (3.10), namely, the dilatation (zi = azj) and the translation (zi = Zi+fl), 
preserve the type A' space (3.8). The former is trivial while the latter is un- 
derstood from the transformation of cr k under the translation: 

k 

a k (z) h+ a k (z + /3) = J £/3 k - l C l (M)a l (z), (3.12) 

1=0 

where Ci are constants depending on M. We now easily see that the latter 
transformation indeed preserves all the elements of (3.8) within the space. 



4 Construction of Quasi-solvable Operators 



In this section, we shall construct the general quasi-solvable operators of (at 
most) second-order leaving the type A' space (3.8) invariant. The set of linearly 
independent first-order differential operators preserving the type A' space is 
given by, 

M Q 

F^k^U^Q^ (k = l,...,M), (4-la) 
F w = oi (V - 1 - £ k(J k^A ■ (4.1b) 

\ k=l k/ 
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In Eq. (4.1a), {m^ is an abbreviation of the set of M non-negative integers 
defined by 

f M 1 

{rrii}k = Ui, . . . ,m M m { G Z> , ^irrii < k>. (4.2) 

^ i=i J 

In the single variable case (M = 1), the set of differential operators (4.1) 
consists of 



(4.3) 



and hence is essentially the same as the st(2) generators. It is as expected 
since the single- variable type A' space is nothing but the sl(2) module (3.6). 
In the case of two variables (M = 2), the set of differential operators (4.1) is 
composed of 

F {01}2 , 2 = cr 2 — , F {mi0}2t2 = a? 1 — (mi = 0,1, 2), (4.4) 

and is essentially the same as the generators of gl(2) kR 3 in Ref. [19]. Actually, 
the two-variable type A' space (3.8) coincides with the g[(2) x M 3 module 
investigated in Ref. [19]. In this sense, the type A' provides a natural extension 
of that Lie-algebraic scheme of two variables to arbitrary number of variables. 

The set of linearly independent second-order differential operators leaving the 
type A' space invariant is as follows: 

M o 2 

F {mi}k+l M^W^ H Q^ (k,l = l,...,M;k>l), ( 4 - 5 ) 

i=l & ^ 

/ M Q \ M _ Q 

FioF^^a^-l-^—J^ar^ (k — 1, . . . , M), (4.6) 
F 10 F,„ = <tJ (^ - 2 - g to t — J (AT - 1 - g /.,-) , (47) 



/ M o \ / M 

^20,00 =(72 (^-2-^^—1(^-1-^^ 

\ /c — 1 / \ / — 1 



d_ 



(4.8) 



where in Eq. (4.6) {mi}k is an abbreviation of the set of M non-negative 
integers defined by 

{^j} fe = jmi,...,m M ?T7.j e Z> , ^irrii — k>. (4.9) 

i=i J 
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Several remarks are in order. First of all, it is important to note that the set 
of the quadratic form of the first-order operators Fs m .\ k ^ in Eq. (4.1a) cannot 
exhaust the second-order operators of the form given by Eq. (4.5) when the 
number of the variables is greater than two (M > 2). In the case of M = 3, 
for instance, the following operator which is an element of Eq. (4.5) 

F <-»W'=^9^? <410) 

cannot be represented by any quadratic combination of the first-order op- 
erators in Eq. (4.1a). Second, even when the number of the variables is two 
(M = 2) where the type A' space (3.8) provides a g[(2) kR 3 module, there exist 
higher-order operators preserving the type A' space which cannot be expressed 
as a polynomial in the first-order operators (4.4), such as Eq. (4.8). In this 
respect, see also Ref. [22]. Third, the operators F w F{ mi } k ^ with Yl^=i^ m % < k 
are represented by linear combinations of the other operators (4.1)-(4.8) as 

/ M \ M o M M 02 

/ M \ M 

\ 1=1 / 1=1 

where {m'^k is obtained from {m^ by m- = m, + while {m'-} k+ i is from 
{rrii}k by m'[ = rrii + 5i : i + 5^. The restriction J2iLi^ m i < ^ ensures that 
E^i ^ < and J2iii irn" < k + I. Thus, F {m ^ ktk and -F{ m »} fc+iifci are in fact 
members of the operators in Eqs. (4.1a) and (4.5), respectively. Hence, the 
operators F 10 F{ rn .y kk with Y^iL\i> m i < k are linearly dependent on the other 
operators (4.1)-(4.8). That is why in Eq. (4.6) we restrict the values of the set 
{rrii} to the one given by Eq. (4.9). Furthermore, the operators of the form 
F{mi}k,k F io are also linearly dependent on the set of operators (4.1)-(4.8) since 
the anti-commutator of F^ m ^ k ^ k and Fi reads 



F(rn. 



{mi} k ,ki -^10 



k,iF w . (4.12) 



Therefore, the most general quasi-solvable operator of (at most) second-order 
which preserves the type A' space (3.8) is given by the linear combination of 
the operators (4.1)-(4.8): 

M M 

T^-M =~X] ^-{mi} k+ i,kl F {mi} k +l,kl - X] ^10,{m}k,kFloF{rn,} k ,k 

k>l {nii} k+l k=l {mi} k 

— ^10, 10-^10-^10 — ^20,00-^20,00 

M 

+ E E B {mi} k ,kF {mi } k ,k + B W F 1Q - c , (4.13) 

k=l {nii} k 
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where the coefficients A, B with indices and Co are real constants and the 
summation over the set {m^ etc. is understood to take all the possible set 
of values {mi,...,mjif} indicated in Eqs. (4.2) and (4.9). In terms of the 
variables a, the operator Tiff ■* is expressed as 

M d 2 
H% ) = " E [M<r)M<Tk<ri + A kl (a' ' 

k,l=l 

M g 



da k d<Ji 

+ £ [B (a)ka k - B k {a)\ - C(a), (4.14) 



k=l 



where A , A kh B , B k , and C are polynomials of several variables given by 
A (a) = A 10AO al + A 20 ,ooO"2, (4.15a) 

M M 

A kl (a) = A {m i } k+l ,kiU' T T i ~ E A 10j{mi}kjk la 1 <j l ]l<j™ i , (4.15b) 

B (<t) =(2Af-k- 3)A (a) - B 10 a u (4.15c) 

M M 

B k {a)={M-k-l) £ ^lo.^H^iIl^r- E S Wfc> *nO, 

(4.15d) 

C(<r) ={N-l){N- 2)A (a) - (Af - 1)Sio<ti + c . (4.15e) 

Among the set of the operators (4.1)-(4.8), ^{m;}^*: and -F{ m4 } fe+; ,M preserve the 
type A' space (3.8) for arbitrary natural number J\f. Therefore, the operator 

~ (A') 

Tt\f is not only quasi-solvable but also solvable if 

Aw,{mi} k ,k = A 10jl0 = A 2 ofio = B 10 = 0. (4.16) 



5 Extraction of Schrodinger Operators 



In the preceding section, we have constructed the most general quasi-solvable 
second-order operator Ttj^ preserving the type A' space (3.8). By applying 
a similarity transformation on Tiff ■* and a change of variables, we may ob- 
tain a family of quasi-solvable operators of a desired form. However, second- 
order linear differential operators of several variables are in general not gauged 
Hamiltonians, that is, they cannot be always transformed to Schrodinger op- 
erators. This fact is one of the most obstacles in constructing quasi-solvable 
quantum many-body systems. Recently in Refs. [11,18,23], it was shown that 
the amount of the difficulty can be significantly reduced by considering the 
underlying symmetry of the invariant space of quasi-solvable operators. In the 
present case, however, there is no full GL(2,M.) invariance, especially no in- 
variance under the special projective transformation, as we have mentioned 
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previously in Section 3. It turns out that the dilatation and translation in- 
variance are insufficient to extract the general form of the type A' gauged 
Hamiltonians, namely, the operators which can be transformed to Schrodinger 
operators from the most general type A' quasi-solvable second-order operators 
(4.14). 

Let us first review the general condition under which a second-order linear 
differential operator of several variables can be cast into a Schrodinger operator 
[23]. Suppose the operator under consideration 7i has the following form in 
variables z: 

M r) 2 M r) 

i,j=l u^uzj i=1 uz t 

Then, it is readily shown that Ji can be cast into a Schrodinger operator by a 
gauge transformation and a change of variables Z{ = z{qi) such that 



c 



W ^ W = ~f:^ + V(lh (5-2) 



2 U a?. 

if and only if the following conditions are satisfied: 

P ij (z) + P ji (z) = (i>j), (5.3) 
{z[f = 2P u (z), (5.4) 
™ _ S t (z) z'l 

where z[ denotes the derivative of zi with respect to The first condition 
(5.3) in general consists of a set of algebraic identities. The second and third 
conditions (5.4)-(5.5) on the other hand are sets of differential equations and 
do not necessarily have a solution. In order that the second condition has a 
solution, the r.h.s. of Eq. (5.4) must depend only on the single variable zf 

(z' i ) 2 = 2P u (z)=2A(z t ), (5.6) 

since we have assumed that the change of variables is determined by a single 
function of a single variable Z{ = z{qi) and thus the l.h.s. of Eq. (5.4) depends 
solely on the variable Furthermore, in order that the third condition (5.5) 
has a solution, the following integrability condition must be fulfilled for all 

i ^ I: 

d dW d dW 1 dSAz) 1 dSAz) 

777s a. = 777s 51 » ( 5 - 7 ) 



dqi dqi % dqi A(zi) dz x A(z t ) dz, t 

where Eq. (5.6) is employed. Therefore, the operator TC can be transformed 



10 



to a Schrodinger operator if it has the form of 

M q2 M q 

n = - £ + E - n*), (5-8) 

i=l aZ « i=l 

with A(zi) and S^z) satisfying Eqs. (5.6) and (5.7). 

Next, we note that in our present case we have made the additional change 
of variables from z to a, Eq. (3.5). From the formula (A. 4), the part of the 
second-order operators in Eq. (4.14) has the following form in terms of z: 



M Q2 M 02 

E [M°)Mwi + A »( ff )]^r^r = - E p «Wor«7. (5 - 9) 

k,l=l UOkUOl ij=1 UZiUZj 



with 

_ , x E^ = i(-l) fc+ ^f- fe #-'[A (a)^^ + A„(a) 



l JV ' TfM /„ „ \ttM 



n m (^i)(^i — 2jn) nn(^)(^j — ^n) 



(5.10) 



To satisfy the condition (5.6), it is evident that the denominator in the r.h.s. 
of Eq. (5.10) for Pu(z) must be completely canceled with a part of the nu- 
merator; otherwise, Pu(z) cannot be a function of the single variable Z{. As a 
consequence, the part of the second-order operators in Eq. (4.14) which satisfy 
the conditions (5.3) and (5.6) must have the following form: 

M a2 M I n \ 02 

-£ [AWHw + Aw(ff fe = -§(5 vf W (M1) 

where constants. Then, the next task is to examine which of dp Cclll DG 

non-zero free parameters. From a dimensional analysis, we easily see that the 
each term in the r.h.s. of Eq. (5.11) for a fixed p must be expressed in terms 
of a as 



Y, z fg^= E E A {mh +l+P - 2 ,ki II °T l q d > ( 5 - 12 ) 

i=l °^ fc,i=i {Tn i } fc+(+p _ 2 i=l CWT fc a<Tj 



with a suitable set of constants AfeL , , ,. On the other hand, we can read 

from Eqs. (4.15a)-(4.15b) that the l.h.s. of Eq. (5.11) contains operators of 
the form 

M Q2 

n^ m< ^-, (5.13) 

only for < Y^iLi i> m i < k+l+2. From this fact and Eq. (5.12), we readily know 
that Eq. (5.11) cannot be satisfied unless a p = for p > 4. From Eq. (4.15b), 
we see that the l.h.s. of Eq. (5.11) contains all the set of operators of the form 
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(5.13) for (0 <)k + I - 2 < J2Zi im i < k + l, and hence Eq. (5.11) can be 
satisfied for all non-zero values of oo, a>i, and a 2 , with suitable values of the 
parameters A{ m .} k+hk i in A k i(cr). However, only the second term in the r.h.s. of 
Eq. (4.15b) produces the operators of the form (5.13) for YhL\ i m i — k + l + 1 
and is insufficient to express the operator in the r.h.s. of Eq. (5.11) for p — 3, 
cf. Section B.5 in Ref. [11]. Hence, Eq. (5.11) cannot be satisfied unless 03 = 
and Aio,{m:i} fe ,fc = 0. A similar observation for Y^iLi = k + l + 2 leads to 04 = 
and Ai ,io = A 2 o,oo = 0. Summarizing the analyses so far, we conclude that 
a p = for all p > 2 and thus the coefficient of the second-order operator of 
the type A' gauged Hamiltonians of the form (5.8) must be 

A(zi) = a 2 zl + a±Zi + a . (5.14) 

Simultaneously, all the following parameters in A (a) and A k i(a) must vanish: 

AlO t {m t } k ,k = -^-10,10 — ^20,00 = 0. (5.15) 



Next, we shall examine the part of the first-order operators. With the aid of 
Eq. (A. 4), the part of the first-order differential operators in the general type 
A' operators (4.14) reads 



M a M o 

£ [B (a)ka k - B k (a)} — = ^ S t (z) — , 

k=l 



da k 



i=i 



dzi 



(5.16) 



with 



Si(z) 



Ef =1 (-l) fc+1 ^ M - fe [B (a)^ fe - B k (a) 



Si(z) 



(5.17) 



Hjf^i) ( z i - zj ) nlf^i) fa - zj ) 

Thus, the derivative of Si(z) with respect to z\{l 7^ i) in our case is calculated 



as 



dSj(z) _ (zj - zi)diSj(z) + Sj{z) 

dz i ~ ziyn^izi- z^ 



(5.18) 



Therefore, the integrability condition (5.7) can be satisfied if and only if the 
numerator of the r.h.s. of Eq. (5.18) has the following form: 

M 

(zi - zi)diSi(z) + Si(z) = fi(zi, zi)f 2 (a)A(zi) J| (^ - zj), (5.19) 

where fi(zi,zi) = fi(zi,Zi) and /2(c) is a function depending solely on the 
elementary symmetric polynomials. 



On the other hand, under the condition (5.15) satisfied, the coefficient of the 
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first-order differential operators in Eq. (4.14) reads 

M 

B (a)ka k - B k (a) = -B 10 ka 1 a k + £ B {mi}kjk ]J <j™\ (5.20) 

{mi} k i=l 

Hence, Si(z) defined by Eq. (5.17) is a polynomial of degree at most M + 1 
in the variables z. From Eq. (5.19), we conclude that the combination of the 
functions f\(zi, zi) f2((r)A(zi) in the r.h.s. of Eq. (5.19) must be a polynomial 
of degree at most 3 in the variables z. Let us first examine the highest- degree 
term. It comes from the first term in the r.h.s. of Eq. (5.20). The corresponding 
term in Si(z) is 

M 

S t (z) = -B 10 a 1 ^(-l) fc+1 ^f-V fc . (5.21) 
k=i 

With the aid of the formulas (A.1)-(A.3), we have 

M 

(zi - zi)diSi(z) + Si(z) = -B 10 Zi(zi - zif Yl ( z i - z j)- ( 5 - 22 ) 

Comparing with Eq. (5.19), we see that f\(zi,zi) = (zj — z\) 2 and fi{<J) = 
const, for the highest- degree term. However, the remaining term — B w Zi is a 
monomial of first-degree. This means that the highest- degree term (together 
with lower-degree terms) cannot be expressed as Eq. (5.19) unless A(zi) is a 
polynomial of first-degree. In other words, the highest- degree term can exist 
if and only if a<i = in Eq. (5.14). Hence, the possible form of the formula 
which includes the highest-degree term must be 

M 

(zi - zi)diSi(z) + Si(z) = -2gA(zi)(zi - zi) 2 5 a2fi Yi ( z i ~ z i)i ( 5 - 23 ) 

where g is a constant. It is evident that the term proportional to a\ in the 
r.h.s. of Eq. (5.23) comes from the highest-degree term given by Eq. (5.22) 
and thus the constant B 10 in Eq. (5.22) is expressed as 

B 10 = 2g ai 5 a2fi . (5.24) 

We will later see that there in fact exists the term in Si(z) which corresponds 
to the term proportional to a in the r.h.s. of Eq.(5.23). Next, we shall ex- 
amine the lower-degree terms, namely, the terms of degree less than M + 1 
in Si(z) which result in the terms of degree less than 3 in the combination 
fi(zi, zi)f2(cr)A(zi). As we have already obtained the terms which can only 
exist in the case of a>i = 0, we can assume here that ai 7^ and thus A(zi) 
is a strictly second-degree polynomial. Hence, only the possible form for the 
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lower-degree terms is fi(zi, zi)f2(cr) = — 
(zi - zi)diSi(z) +Si{z) = - 



2c, c is a constant, and we have 



M 

-2cA( Zi ) J] (zi-Zj). 



(5.25) 



Substituting the higher-degree term (5.23) together with the lower-degree term 
(5.25) in Eq. (5.18), we eventually obtain 

^^ = -2gA(z t )5 a2fi -2c T ^- (i^l). (5.26) 

OZi yZi Zi) 

This set of differential equations can be easily integrated as 

Si(z) = -Q( Zi ) - 2ga 1 A(z i )S a2>0 - 2c ]T (5.27) 

where Q{z,j) is a function of a single variable. This term should come from 
at most Mth-degree terms in Si(z) which cancel the denominator of the last 
term in Eq. (5.17) so that it depends solely on the single variable Zj. The 
denominator is of degree M — 1 and thus Q{zj) must be a polynomial of at 
most first-degree: 

Q( Zi ) = b lZi + bo. (5.28) 

Next, we shall check whether Eq. (5.16) actually holds for Eqs. (5.20) and 
(5.27). From the formulas (A. 5), the r.h.s. of Eq. (5.16) with Si(z) given by 
Eq. (5.27) is expressed as 

m q m r 

Y / S l (z)— = J2\(M-k + l)(M-k + 2)ca a k _ 2 

i=l aZ i k=l ^ 

- (M - k + 1) [&o + (M - fc)cai] <r fc _i - k [&i + (2M - fc - l)ca 2 ] cr fc 

1 «9 

- 2(M - fc + 1)5-00^2,00-1^-1 - 2a 1 g5 a2:0 ka 1 a k j— . (5.29) 

We now easily see that all the terms in the braces in the r.h.s. of the above 
equation are contained in Eq. (5.20) for each fixed k. Hence, all the term in 
Eq. (5.27) can indeed exist. Substituting Eq. (5.27) in Eq. (5.5) and using 
Eq. (5.6), we have 

dW Q( Zi ) A'(zi) c ™ 1 . 

~« = "t77 \ + T77 V — 9< T i°a 2 ,o ~ c • (5-30) 

dzi 2A(zi) 4:A{Zi) ^ Zi - Zj 
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Again, we can easily integrate the above set of differential equations to obtain 



M . 

w = - E J dz 



Qjzj) 
2A( Zi 

— c y^ln \zj — z. ' 



1 



M 



i=l ' 
M 



4 1 



i=i 



7 b 



(5.31) 



where we have omitted the integral constant. 



Finally from Eqs. (5.8) and (5.27), we find that the gauged HamiltoniansiJy^ ^ 
which preserve the type A' space (3.8) must have the following form: 



Er(A') 



M 



H 



d 



M 



d 



M 



AT 



A( Zi ) a 



dzi 



%i Zj dz{ 



1=1 ^"l 1=1 

where A{zA and B^z) are given by Eq. (5.14) and by 

Bi(z) = 2ga 1 A(z i )5 a2fi + Q(zi), 



C(a), (5.32) 



(5.33) 



respectively. From Eqs. (4.15e), (5.15), and (5.24), the function C in Eq. (5.32) 
is calculated as 



C(a) = -2(JV - l)#ai(7i<5 a2) o + c . 



(5.34) 



It is easily shown that the gauged Hamiltonian (5.32) can be actually cast into 
a Schrodinger operator by a gauge transformation 



M 



d 2 



1 



M 



2 g% 2 + 2^ 



i=l 



9<H , 



d 2 W 
~W 



(5.35) 



if the gauge potential W is chosen as Eq. (5.31) and the function z(q) which 
determines the change of variables satisfies 



z\qf = 2A(z(qj) = 2(a 2 z(qf + a lZ (q) + a ). 



(5.36) 



In this case, Eqs. (4.16), (5.15), and (5.24) tell us that the gauged Hamiltonian 
(5.32) is not only quasi-solvable but also solvable if and only if 



B w = gai5 a2fi = 0. 



(5.37) 



It is apparent from the construction that the Hamiltonian (5.35) preserves the 
space Vjf.M defined by 

V#£ = e- w V#i- (5.38) 
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Hence, the type A' Hamiltonian ^ can be (locally) diagonalized in the 
finite dimensional space (5.38). We shall thus call the space (5.38) the solvable 
sector of Hj^ \ 



6 Classification of the Models 



We shall now explicitly compute the concrete form of the type A' quan- 
tum Hamiltonians. From Eqs. (5.14), (5.30), and (5.34), the potential term 
in Eq. (5.35) is explicitly calculated in terms of z as 



M 



y = J2 



+ 9 



M 



2Q(z t ) - A'(zi)\ [2Q( Zi ) - 3A'( Zi ) 

M 



M 



a^z) £ Q( Zi ) + (ga^zf + l) £ A( Zi ) + a ± (M, A/") £ 

i<j \ z i z i) 



(6.1) 



where, and in what follows, V denotes an arbitrary constant, and the coupling 
constant a 1 (M,A/') is given by, 



ai (M,A0 = [2(Af- 1) + M(M- l)c]ai. (6.2) 

From Eq. (5.36), the change of variable is determined by the following integral: 

dz 



±(q -qo)= — f 



2(a 2 z 2 + a 2 z + a x ) 



(6.3) 



In contrast to the type A case where the systems are constructed from the 
sl(M + 1) generators, our present models do not have full GL(2, R) invariance 
as has been mentioned previously in Section 3. However, the remaining invari- 
ance under the dilatation and translation enables us to classify the type A' 
models. Indeed, it is readily shown that A(z) can be cast into one of the canon- 
ical forms listed in Table 1 by a combination of the dilatation and translation. 



Furthermore, we note that from Eq. (5.36), a rescaling of the coefficients Oj, 
bi, Co by an overall constant factor v has the following effect on the change of 
variable z(q): 

z(q ; vai, ubi, vcq) = z(y/uq; a u b h c ). (6.4) 
From this equation and Eqs. (5.14), (5.28), (5.31), and (6.1), we easily obtain 
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Table 1 

Canonical forms for the polynomial A(z), (5.14), where v is a positive real constant. 



Case 


Canonical Form 


I 


1/2 


II 


2z 


III 


±2uz 2 


IV 


±2v{z 2 - 1) 


IV 


±2v{z 2 + 1) 



the identities 

W(q;uai,ubi,uco) = W{^i> q ; a u b h c ), (6.5a) 
V(q ; va h vb h uc ) = v V{^/vq ; Oj, h n c ). (6.5b) 

We shall therefore set v — 1 in the canonical forms in Cases II-lV'), the models 
corresponding to an arbitrary value of v following easily from Eqs. (6.4) and 
(6.5). It should also be obvious from Eq. (6.3) that the change of variable 
z(q), and hence the potential V determining each model, are defined up to 
the transformation q i— > ±(q — q ), where q G K. is a constant. The solvability 
condition (5.37) implies that except for the model with g ^ corresponding 
to Case II in Table 1 all the obtainable models are not only quasi-solvable but 
also solvable. 

As we will see below, the potentials in Case I, III, and IV have singularities 
at qi = qj for all i ^ j in the subspace {— Q < qi < Q} G 1R M where Q is a 
submultiple of the real period of the potential or Q = oo when the potential 
is non-periodic. Similarly, the potentials in Case II and IV are singular at 
qi = and q t = ±qj for all i ^ j in the same subspace of M M . Hence, the 
Hamiltonians are naturally defined on 

< q M < ■ ■ ■ < qi < ^. (6.6) 

Normalizability of the solvable sector (5.38) on the space (6.6) depends on 
the behavior of the gauge potential W(q) in the each case. The finiteness 
of the L 2 norm of the two-body wave function in the solvable sector Vj^.m m 
general leads c > —1/2, where c denotes the coupling constant of the two-body 
interaction appeared in the last line of Eq. (6.1). 



17 



6.1 Case I: A(z) = 1/2 



Change of variable: z(q) = q. 



Potential: 



(M \( M \ M 1 M 

V(q) = 9 1^-9 + bij (E Qij + M bog E 9< + ^ E( 6 ift + 6 °) 



+ c(c-l)i: ¥ ^ + ^ (6.7) 



Gauge potential: 



i<j 



M \ 2 l M M M 



w(?) = -f(Eft) -yE^-^Eft-cE^lft-^l- ( 6 - 8 ) 

\i=l / i=l i=l i<j 

When (7 = 0, this case corresponds to the rational A M -i type Calogero- 
Sutherland model [5] and is identical with the solvable model corresponding 
to Case I of the type A models, Eqs. (7.8)-(7.9) with b 2 = in Ref. [11]. 
Hence, the above model provides an example of deformed Calogero-Sutherland 
models which preserve quantum solvability. The parameter 6 corresponds to 
the translational degree of freedom and is irrelevant. In fact, if we employ 
the translational freedom mentioned earlier below Eqs. (6.5) and apply the 
translation qi i— > qi — b / {b\ +Mg) in Eqs. (6.7)-(6.8), the potential and gauge 
potential become 



V(q) 



(M \ ( M \ 2 b 2 M M 1 



<3 

(6.9) 

( M \ 2 l M M 

w(«) = -f Eft -^E^-cE ln lft-^i (6-io) 

z \i=l / z 1=1 i<j 

and have no dependence on 6 any more. It is now obvious that the model is 
exactly the rational Am-i type Calogero-Sutherland model with the center- 
of-mass coordinate subjected to the harmonic oscillator potential. Since we 
can easily separate the center-of-mass coordinate from the others, quantum 
solvability of the above model is readily understood. 

The one-body part of the potential has no singularities and hence a natural 
choice is Q = oo. In this choice, the form of the gauge potential (6.8) tells us 
that the solvable sector (5.38) is square integrable on the space (6.6) as long 
as g < 0, bi < 0, and c > —1/2. Hence, the model (6.7) is exactly solvable on 
Eq. (6.6) in these parameter regions. 
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6.2 Case II: A(z) = 2z 



Change of variable: z(q) = q 2 



Potential: 



M 



V(q) = (j + 2(2AT - l)g + Mb g + 2M(M - l)cgj £ q\ 



M 



M 



+ 2, 2 (£g 2 ) + M(£ft 2 J + 
1 1 



\i=l / 
M 

+ c(c-l)E 

i<j 



(ft - ft) 2 (ft + ft) 2 . 



(fe _l )(6o _ 3) M i 
8 ^ft 2 



+ ^0- 



(6.11) 



Gauge potential: 



_ / M \ 2 L M 



M 



7> - 1 M 

^Elnlftl-cEHft 2 -^. (6.12) 



This case corresponds to the rational -BCm type Calogero-Sutherland model 
[24] when g = 0. Furthermore, the above model is quasi-solvable when g ^ 0. 
To the best of our knowledge, it is new and provides another quasi-solvable 
deformation of Calogero-Sutherland models which is different from the In- 
ozemtsev type classified in Ref. [11]. In this respect, we should refer to sim- 
ilar quasi-solvable models in the literature, namely, Eq. (23) in Ref. [7] and 
Eq. (3.7) in Ref. [25]. These models are deformations of the rational Am-\ 
type Calogero-Sutherland system and thus their two-body interaction terms 
are different from our BCm type. Furthermore, the solvable sectors of these 
models are spanned by monomials of a single variable while that of our present 
model is by monomials of M variables (3.8). 

The one-body part of the potential is only singular at qi = and hence a 
natural choice is f2 = oo. In this choice, we see from Eq. (6.12) that the 
solvable sector (5.38) is square integrable on the space (6.6) as long as g < 
and c > —1/2. Hence, the model (6.11) is quasi-exactly solvable on the space 
(6.6) in these parameter regions. When M — 1, the above model becomes 



V{q) = \q\Agq 2 + + (4Af + b - 2)gq 2 + (& ° ^ 3) + V , (6.13) 

and thus exactly reduces to the well-known one-body quasi-solvable sextic 
anharmonic oscillator, classified in Case II of the type A models, Eqs. (7.11)- 
(7.12) with b 2 = Ag and c x = (b - l)/2 in Ref. [11]. 
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6.3 Case Ilia: A{z) = 2z 2 



Change of variable: z(q)=e 2q . 



Potential: 

h (h A\ M h 2 M M 1 

(6.14) 

Gauge potential: 

= T E e" 2 * - (| - 1 + (M - l)c) £ 
4 i=i v z y i=i 

M 

— c^ln|sinh(gi — qj) . (6.15) 

i<j 

This model is the hyperbolic Am-i Calogero-Sutherland model [6] in the 
external Morse potential and completely the same as the model corresponding 
to Case III of the type A models, Eqs. (7.15)-(7.16) with b 2 = in Ref. [11]. 



6.4 Case Illb: A(z) = -2z 2 



The formulas of the potential and gauge potential for this case can be easily 
reduced from Eqs. (6.14)-(6.15) using Eqs. (6.5) with v = —1. The change of 
variable is z(q) = e 2iq . 



6.5 Case IVa: A(z) = 2(z 2 - 1) 



Change of variable: z(q) = cosh2g. 



Potential: 



V(q) 



(6 1 -6 -2)(6 1 -6q-6) 
8 

M 

+ c(c-l)E 



M 



\ sinh 2q, 



+ 



6o(6i-4) 



M 



+ 



1 



sinh (qi - q 5 ) sinh + qj) _ 



^ sinh 2 qi 
+ V . (6.16) 
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Gauge potential: 



6i-2 



M 



b 



M 



W(q) = — ^2 m I smn 2ft | — -^ln| tanh q. 



i=l 



i=i 



M 



c^ln sinh(ft - ft) sinh(ft + ft) 



(6.17) 



i<j 



This model is the hyperbolic BCm Calogero-Sutherland model [24] and com- 
pletely the same as the model corresponding to Case IV of the type A models, 
Eqs. (7.22)-(7.23) with b 2 = in Ref. [11]. 



6.6 Case IVb: A(z) = -2(z 2 - 1) 



The formulas of the potential and gauge potential for this case can be easily 
reduced from Eqs. (6.16)-(6.17) using Eqs. (6.5) with v = —1. The change of 
variable is z(q) = cos2g. 



6. 7 Case IV a: A(z) = 2(z 2 + 1) 



Change of variable: z = sinh2ft 



Potential: 



v{q) _ bl-(b 1 -2)(b 1 -6) ^ : 



M 



i<j 



+ 'v- - 1) £ :~:2, — - —2 



1 | 6o(6i-4) ^ sinh 2ft 

cosh 2 2qi 4 i=1 cosh 2 2ft 

1 1 



sinh (qi - ft) cosh (ft + ft)_ 



+ V . 



(6.18) 



Gauge potential: 



6i-2 



M 



X] In I cosh 2 ft I - T E§ d 2 1i 



i=i 



i=i 



M 



sinh(ft - ft) cosh (ft + ft) 



(6.19) 



where gdg = arctan(sinh g) is the Gudermann function. This model is another 
hyperbolic BCm Calogero-Sutherland model and completely the same as the 
model corresponding to Case IV of the type A models, Eqs. (7.26)-(7.27) with 
6 2 = in Ref. [11]. 
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6.8 Case IV b: A(z) = 



-2(z 2 + 1) 



The formulas of the potential and gauge potential for this case can be easily 
reduced from Eqs. (6.18)-(6.19) using Eqs. (6.5) with v = —1. The change of 
variable is z(q) = isin2g. 



7 A New Type C Generalization Based on the Type A' Scheme 

In the preceding sections, we have investigated (at most) second-order linear 
differential operators preserving the type A' space (3.8). This new generaliza- 
tion of the single-variable type A space to several variables suggests a new 
multivariate generalization of the single-variable type C space. The type C 
monomial space of a single variable z is defined by [17] 

V& = V^®^V^, (7.1) 

where (i = 1,2) is a type A monomial space of dimension Mi defined by 
Eq. (3.6), Mi and M2 are positive integers satisfying Mi > M2, and A is a real 
number with the restriction 

XeR\{-M 2 ,-M 2 + i,...,Mi}, (7.2) 

and with A 7^ —M2 — 1, Mi + 1 if Mi — 1 or M2 = 1. In the previous paper [18], 
we generalized the type C space of a single variable (7.1) to several variables 
as follows: 

V M%N%M — v Hr,M v M2\Mi [i- 6 ) 

where Vjfi} M (i — 1, 2) is a multivariate type A space defined by Eq. (3.7). The 
above generalization scheme based on the type A space (3.7) now strongly sug- 
gests another generalization based on the type A' space (3.8) as the following: 

Indeed, the latter space also reduces to the single- variable type C space (7.1) 
when M — 1. Hence, Eq. (7.4) would provide a new multivariate generalization 
of the type C space and we hereafter call the space (7.4) type C. In this section, 
we shall investigate (at most) second-order linear differential operators which 
preserve the type C space (7.4). 

First of all, the set of linearly independent first-order differential operators 
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preserving the type C space is given by, 



M d 

w*=n<v (*=i,...,m-i), (7.5) 

i=l acr fe 

^jmm = ff*f7 — • (7.6) 

Similarly, the set of the linearly independent second-order differential opera- 
tors preserving the type C space is as follows: 

^=11^^ (fc,r=l,...,M-l;fc>Z), (7.7) 

i? w 5^ = J] «— «- (* = 1,...,M-1), (7.8) 

m did \ 
F {mi]M ,M{E MM - A) = J " A j ' (7 ' 9) 

/ M <9\/ M <9 \ 

f 10 , 00 s - 1 - - 1 -X>^} ("o) 

Therefore, the most general quasi-solvable operator of (at most) second-order 
which preserves the type C space (7.4) is given by the linear combination of 
all the operators (7.5)-(7.10): 

M-l M-l 
ft(C) = _ ^ ^ ^{mi}fc +r ,kl F {mi}- k+T ,M ~ J2 A {rn l }- k AMMF {mi} -- k E MM 
k>l {mi}k + i k=l {m,i} % 

— E A{ mi y M! M,MMF{ mi y M! M(EMM — A) — A 10i00 F 10i oo 
{"ii} M 
M-l 

+ E B {m z }- k ,kF{ mi }- k ,k + BmmEmm - c , (7-11) 

k=l 

where again the coefficients A B with indices and Co are real constants and 
the summation over the set {rrii} k etc. is understood to take all the possible 
set of values {mi, . . . , m M } indicated in Eq. (4.2). In terms of the variables a, 
the operator 7v c ^ is expressed as 

-ff[C>) = _ Ui ,oofcZ<7i<7 fc <7 i + A kl {a)] 

k% 1 ^o k da x 

+ E [^io,oo(AT + MX - k - 2)ka l a k + B k (a)] — - C(a), (7.12) 



k=i 
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where Am, Bk, and C are polynomials of several variables given by 



M 

A M (a)= £ ^WwIlO (M = l,---,M-l;fc>Z), (7.13a) 

M 

A m {a) = A {mi}k , k ,MM II a T i(J M (k = l,...,M), (7.13b) 

{m.;} fe i=l 
M 

B S ((7) = - E IK' (* = 1, . . . , M - 1), (7.13c) 

M 

Bm{o) = —(A — 1) A{ m .} mMtMM Y[ — Bmm&m, (7.13d) 

{mi} M i=1 

C{a) = (M - 1)(JV 2 + MA - l)A 10 ,oo^i + c . (7.13e) 

Except for the operator F 10i00 , all the operators in Eqs. (7.5)-(7.9) leave the 
type C space (7.4) invariant for arbitrary natural numbers M\ and J\f 2 . Hence, 
the operator W c '^ is not only quasi-solvable but also solvable if 

^io,oo = 0. (7.14) 



As in the case of regular type C in Ref. [18], all the operators in Eqs. (7.5)- 
(7.10) and hence the most general type C operator (7.12) preserve separately 
both the subspaces of Vj^ jj t m m Eq- (7.4), the fact originally comes from the 
restriction (7.2). In other words, the second-order operators a^ k ~^ x iiS c '^ a^~^ x 
(k = 1,2) leave the type A' space vffj invariant, respectively. As a conse- 
quence, the most general type C gauged Hamiltonian H^ c "> must satisfy the 
following condition: 

In the above, each of the operators H\fJ (k = 1,2) is a type A' gauged 
Hamiltonian and thus has the form (5.32) with Eqs. (5.14), (5.28), (5.33), and 
(5.34). Tracing a completely similar way to Section 5 in Ref. [18] and noting 
the fact that the first-order operator Fiq defined by Eq. (4.1b) is missing in 
the set of the type C operators (7.5)-(7.10), we find that the most general 
type C gauged Hamiltonian satisfying the condition (7.15) has the following 
form: 

M a2 M q M a( \ q 
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where 

A(zi) = a 2 Zi + aiZi, (7-17) 
B(z l )=b 1 z l -(X-l)a 1 , (7.18) 

and cij, bi, cq, and c are constants. Comparing the above results with the most 
general type C gauged Hamiltonian, Eqs. (5.17)-(5.20) in Ref. [18], we see 
that the type C gauged Hamiltonian (7.16) is a special case of the type C 
gauged Hamiltonian with a 3 = 0. Hence, all the quantum mechanical models 
of type C are included in the ones fully classified in Ref. [18]. 



8 Discussion and Summary 

In this article, we have made a new generalization of the type A monomial 
space of a single variable to several variables and have constructed the most 
general (at most) second-order quasi-solvable operator which preserves the new 
linear space called type A'. Examining the condition under which the type A' 
second-order operators can be transformed to Schrodinger operators, we have 
extracted the most general type A' gauged Hamiltonian. Then, we have com- 
pletely classified the type A' quantum Hamiltonians. We have also investigated 
a new type C generalization called type C based on the type A' space. Com- 
bining the results obtained in this article with the ones in Refs. [11,18], we can 
summarize the classification of the type A' quasi-solvable quantum many-body 
systems as shown in Table 2. 

Table 2 

Classification of the type A' (quasi-) solvable quantum many-body systems. 



Model 


Type 


Solvable 


Rational A Calogero-Sutherland 


A', A 


o 


+ quadratic M-body interaction 


A' 


o 


Rational BC Calogero-Sutherland 


C, C 


o 


+ sextic M-body interaction 


A' 


X 


Hyp. (Trig.) A Calogero-Sutherland 


c, c 


o 


+ external Morse potential 


A', A 


o 


Hyp. (Trig.) BC Calogero-Sutherland 


c, c 


o 



The meaning of Table 2 is as follows. For example, the rational BC Calogero- 
Sutherland model in the third row belongs to type C and C, and it is not only 
quasi-solvable but also solvable. If the sextic M-body interaction is added to 
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this model (the fourth row), it belongs to type A' and it is only quasi-solvable 
but not solvable. The interpretation for other models would be straightforward 
in a similar way. 

One may be curious about the fact that some of the models belong to both 
type A and A', or to both type C and C. Let us first consider the former 
cases. The fact that a model belongs to both type A and A' means that it 
preserves type A and A' spaces simultaneously. In this respect, we note that 
for fixed values of M(> 1) and Af the type A' space (3.8) is a subspace of the 
type A space (3.7): 

V#2r C V$ M . (8.1) 

However, it does not mean that an operator which preserves a type A space 
always preserves a type A' space too nor mean vice versa. For instance, all the 
operators of the form 

d 

Eij = CTl da~ ( i >j = 1 >---> M '> i >j)> ( 8 - 2 ) 

preserve type A spaces but do not any type A' spaces, while the operator 
F w defined by Eq. (4.1b) preserves the type A' space with the same Af as in 
F 10 but does not any type A spaces. From the set of the operators (4.1) and 
(4.5)-(4.8) which preserve the type A' space and the set of operators (3.10) 
and (3.12) in Ref. [11] which preserve the type A space 4 , we find that the set 
of linearly independent differential operators of (at most) second-order which 
preserve both the type A and A' spaces simultaneously is the following: 

h a i {l - 3) ' (8 - 3a) 

Q2 ^2 ^2 

«— a~> ai T) — a — (i<k + l), ^a j7 —— (i + j<k + l), (8.3b) 
da k doi da k dai oa k dai 

-^h^ {,+j - k) - (8 - 3c) 




4 We would like to note that there are typos in Ref. [11]; the following operators 
are missing in Eqs. (3.12): 

EoiEoj, 



and Eq. (3.12c) should be 



doidoj 



( M d \ d 

^-1-|>— y 3 — = E t0 E jk . 



26 



The most general quasi-solvable operator of (at most) second-order 7Y^' A '- ) is 
obviously obtained by the linear combination of all the above operators. In 
particular, the most general gauged Hamiltonian H^ A,A ) leaving the type A 
and A' spaces simultaneously invariant must be of the form of both the type 
A gauged Hamiltonian, Eqs. (5.12)-(5.14) in Ref. [11], and the type A' gauged 
Hamiltonian, Eqs. (5.32)-(5.34) with (5.14) and (5.28). From the observation, 
we see that the type A' gauged Hamiltonian (5.32) belongs to type A too if 
and only if g5 a2j o = 0. We note that if it is the case, the system is always 
solvable since the solvability condition (5.37) is automatically satisfied in the 
case. In fact, all the models without M-body interactions in Table 2 satisfy 
the condition g5 a2t0 = 0, thus belong to both type A and A', and are not only 
quasi-solvable but also solvable. One of the interesting consequences is that 
all the models without M-body interactions in Table 2 preserve the following 
infinite flag of finite dimensional linear spaces: 

U U (8.4) 



V {A>) V (A>) y (A') 



where vffi M are gauge-transformed type A spaces 



with the same gauge potential W as in Eq. (5.38). 

The situation in the case of simultaneous type C and C is completely analo- 
gous to in the case of simultaneous type A and A' discussed just above. The 
models which belong to both type C and C preserve the type C and C spaces 
simultaneously. From Eqs. (7.3), (7.4), and (8.1), we easily see that for fixed 
values of M(> 1), Af\, and J\f 2 the type C space is a subspace of the type C 
space: 



The set of linearly independent differential operators of (at most) second-order 
which preserve both the type C and C spaces simultaneously is the following: 
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d 2 

da- k do x 




The most general quasi-solvable operator of (at most) second-order 7i^ c ' c ^ is 
again obviously obtained by the linear combination of all the above operators 
(8.7). It should be remarked that the operator H^ c ' c '^ is always solvable since 
all the operators in Eqs. (8.7) preserve the following infinite flag of finite 
dimensional linear spaces: 

KA/i,l;M <- v Mi,2;M <- <- V Mx^M <- 

U U U (8.8) 

for all M\ = 1, 2, 3, ... . In particular, the type C gauged Hamiltonian (7.16) 
always belongs to type C, as has been mentioned at the end of Section 7, thus 

always preserves the infinite flag of the spaces (8.8) for all M\ = 1,2,3, 

Table 2 indicates that the rational, hyperbolic, trigonometric BC type, and 
hyperbolic, trigonometric A type Calogero-Sutherland systems have this in- 
triguing property. 



Furthermore, it is also interesting to note that the most general operator of 
(at most) second-order 7i^ A ' A ) which preserves both the type A and A' spaces 
for a given M = n always preserves the type A' spaces for all M = 1, 2, 3, . . . 
but does not the type A spaces for any M ^ n as far as the operator given 
by Eq. (8.3c) is included in H^ A '\ As a consequence, the operator H^ A,A '^ 
is always solvable but the infinite flag of finite spaces preserved by it has the 
different structure from Eq. (8.4) as the following: 

y n;M 

U (8.9) 

L V n-1\M L v n\M L y n+l;M L 

Only when the operator (8.3c) does not exist in W A,A '\ the infinite flag of 
finite spaces preserved by it has the same structure as Eq. (8.4). The discussion 
we have made so far surely shows that linear spaces of monomial type which 



V (A ° C ■■■ 

V 1;M <- 
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can be preserved by quantum many-body systems have far richer structure 
than those preserved by one-body systems. The structure we have revealed in 
this article would be just the tip of the iceberg. 

Finally, we would like to recall the fact that the most general type C quasi- 
solvable operator (7.12) preserves separately the subspaces in the type C space 
(7.4) due to the restriction (7.2). It indicates the existence of irregular type 
C operators which do not preserve them separately when the restriction (7.2) 
is omitted, as in the case of type C [18]. These issues on irregular type C 
together with irregular type C would be reported elsewhere. 
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A Formulas 



In this appendix, we summarize some useful formulas on the elementary sym- 
metric polynomials. 

M M 

J2(-l) k z^- k a k = J2(-^M-k = (t=l,...,M), (A.l) 

fc=0 fc=0 
M M 

J2(-l) k+1 kz™ = J] (* - *j) (i = l,...,M), (A.2) 

k=l j(jH) 



dcr k 



k 



£(-*irV*-J (i,fc = l,...,M), (A.3) 



1=1 



F) M ( _-\\k+l y M-k a 

' -E ^M ~ (k = l,...,M). (A.4) 



da k friUf^i)(zi - Zj) dzi 
The first formula (A.l) is readily derived from the following identity: 

M M 

n<* - *) = E(-i)**" _ v 

i=l k=0 

The second one (A.2) is easily proved inductively. The third one (A.3) is 
derived from the repeated application of the following formula: 



da k da k . 

— Cfc-l — z, 



i 



OZi OZj 
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The derivation of the fourth one (A. 4) is as follows. Taking the differential of 
Eq. (A.l), we have 

M ( — l) k+1 z M ~ k 

fc=i L(=o (-l)W-fK 01 

From Eq. (A.l) the denominator in the r.h.s. of the above equation reads, 

Af-l M 

e (-i)'(m - ^r 1 - 1 ^ = E(-i) ,+1 ^" ,_1 ^- 

Hence, applying Eq. (A. 2) we obtain the formula (A. 4). 

For the derivation of the formulas below, see Appendix B in Ref. [11]. 

M Q M Q 

U^rU M - k+1H ^' (A - 5a) 

M Q M Q 

E z iTT = E kcr k^— , (A.5b) 

M 1 q M g 

2 £ = - E( M - * + !)( M - k + 2)^-2^-, (A.5c) 

2E T^VgT = -Y,(M-k)(M-k + iW-ijrr, (A.5d) 

t# Z i ~ Z j ° Z i k =l a<7 k 

M 2 X M a 

%^^r-E k(2M - k - 1)a ^ (A ' 5e) 
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